Black hole pair creation and the stability of flat space 
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We extend the Gross-Perry- Yaffe approach of hot flat space instability to Minkowski space. This 
is done by a saddle point approximation of the partition function in a Schwarzschild wormhole 
background which is coincident with an eternal black hole. The appearance of an instability in the 
whole manifold is here interpreted as a black hole pair creation. 



[f^~. I. INTRODUCTION 

o : 

The issue of stability with respect to quantum fluctuations in the path integral approach to quantum gravity is a 
fundamental problem as important as quantum gravity itself. Due to the attractive nature of the gravitational field, 
C^H' we expect to find a lot of unstable physical configurations. Therefore one might worry about the stability of the 
ground state of quantum gravity. One generally accepts to regard Minkowski space as the ground state, or vacuum, 
of quantum gravity. Classical small perturbations about this vacuum are stable. This is guaranteed by Positivity 

' Theorems of Schoen and Yau0,Qand WittenQ- However one might find that flat space is quantum mechanically 
unstable. This problem was investigated for the first time by Gross, Perry and Yaffe (GPY)Q in the Euclidean 
path integral context. They discovered that hot flat space is unstable and the loss of stability was interpreted as a 
spontaneous nucleation of a black hole with an inverse temperature (3 = 8irMG, where M is the black hole mass and 

■ G is the Newton constant. In the semiclassical approximation to the Euclidean path integral 

f 

Z= Vg^expi-Iglg^)) , (1) 



O 

o : 

the decay probability per unit volume and time T is defined as 

o" 



X 



r = Acxp (-I cl ) = exp {-Ig [g^]) J Vh^ exp [h^ , (2) 



cr 

■ ■ 

?H . where the gravitational field g^ v has been separated into a background field g M „ and a quantum fluctuation h^ u 



Qiiv — diiv "b h^ v , (3) 



A is the prefactor coming from the saddle point evaluation of Z and I c i is the classical part of the action. If a 
\ single negative eigenvalue appears in the prefactor A, it means that the related bounce shifts the energy of the false 
ground state 5]. In the Schwarzschild case GPY discovered one single negative mode \ neg = — .19/ (MG) 2 . Allen 
reconsidered the effect of finite boundaries on the appearance of a negative mode showing that, if the box enclosing 
the black hole is little enough the instability disappears^. Since the pioneering paper of Gross, Perry and Yaffe: 
many other spacetime configurations have been investigated in different contexts. In particular, the de Sitter case 
involving a positive cosmological constant has been examined in the inflationary context by Bousso and Hawking|l8|. 
In the same context of the de Sitter background, Ginsparg and Perry 7], Young 8] and more recently Volkov and 
Wipf 9] have computed the partition function to one loop. On the other hand, the Anti-de Sitter case involving 
a negative cosmological constant has been discussed by Hawking and Page^ an d subsequently by Prestidgepd|. 
The parallel problem in a higher dimensional context of the GPY-instability was firstly taken under examination by 
Witten^^l, who examined the stability of Kaluza-Klein theories. After a decade Gregory and Laflamme reconsidered 
the problem of the gravitational stability in the context of branes[l3j. This opened an interest in string theory where 
the stability of black branes was widely examined by different authors^ . It is interesting to note that the only case of 
instability involving not a nucleation of a single black hole but a pair, is the de Sitter case. Moreover a deep difference 
appears between this case and the Anti-de Sitter or Schwarzschild case: temperatures before and after pair creation 
are different, at least in the neutral de Sitter black hole pair creation. Motivated by this result, we would like to 
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reconsider the stability of flat space without involving temperatures as in the GPY analysis, but looking at a neutral 
black hole pair creation mediated by a wormhole of the Schwarzschild type. Recall that in the GPY instability is 
not possible to reach a vanishing temperature, because T = 1/0 = 1/ (8nMG). Therefore to investigate the stability 
at zero temperature, we need to consider a Schwinger-like process of black hole pair creation 1 . The Schwinger-likc 
process has been investigated in a series of papers in a variational approach based on a Hamiltonian formulation of 
the Einstein gravity In this paper, we would like to apply the GPY analysis to such a problem. The rest of 
the paper is structured as follows, in section [H] we write the gravitational action for the Schwarzschild wormhole, in 
section IiTH we follow the GPY procedure of separating TT modes, in section llVl analyze the eigenvalue equation. We 
summarize and conclude in section Units in which K = c = k = 1 are used throughout the paper. 



II. GRAVITATIONAL ACTION WITH BOUNDARIES FOR THE SCHWARZSCHILD WORMHOLE 



Before considering the perturbations on the Schwarzschild wormhole, we describe some properties of the 
Schwarzschild-Kruskal line clement 



ds 2 = 



32 (MGy 



■ exp 



( 2MGJ 



dU dV + r 2 dn 2 



(4) 



M represents the Arnowitt-Deser-Misner (ADM) mass of the wormhole dfl 2 is the line element of the unit two- 
sphere, and the radial coordinate is regarded as a function of the ingoing and outgoing null coordinates (U, V) as 



uv 



2MGJ 1 V 2MG/ 

while the time coordinate regarded as a function of the same coordinates (U, V) is defined as 

-V 



t = 2MGln 



U 



(5) 



(6) 



The Schwarzschild-Kruskal spacetime is the union of four regions (wedges) R+ , i?_ , T + , and T_ . The regions R+ and 

R- are asymptotically flat . In the R + region, the Kruskal coordinates (U, V) satisfy U < 0, V > 0, while in 

U > 0, V < 0. An important property of the line element (@J) is its invariance with respect to the discrete symmetries 



I : U 



-U: V 



V. 



T : U 



-V: V 



-U 



(7) 



This means that we expect the physics be symmetric with respect to the bifurcation surface So defined as the 
intersection of H + and H~ , which are the future (H + ) and past (H~) horizons, respectively. In Fig.QJ, we show 
the corresponding Penrose diagram We denote with S t the t = constant hypersurface. Any of such hypersurfaces is 

r=0 




r=0 

FIG. 1: This is a Penrose diagram for the Schwarzschild wormhole. 



1 A recent detailed historical overview about the black hole pair creation process can be found in Ref. ll5l 
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invariant under T- reflections and has Einstein- Rosen bridge topology R 1 x S 2 , represented in Fig.©, whose intrinsic 
geometry and time derivatives are chosen to satisfy the gravitational constraint equations We denote by E± the parts 




FIG. 2: Wormhole representation for the Schwarzschild metric. The bifurcation surface So divides the hypersurfaces E+ and 
E_. 



of E lying in R± respectively. The line element (@J restricted to the surface E reads 



dV 



dy 2 +r 2 (y)dn 2 , 



(8) 



where the quantity y, defined by 



dy = ± 



dr 



y/1 - 2MG/r 



(9) 



represents the proper geodesic distance from the "throat" of the bridge located at r = r(y = 0) = r + = 2MG. We 
choose y to be positive in E + and negative in E_. The set (t, y, 9, <p) can be used as canonical coordinates everywhere 
in R± outside the bifurcation surface Sq. These coordinates are right-oriented in i? + and left-oriented in The 
configuration we wish to examine is described by boundary surfaces located in different regions with two different 
slices E' = E t < and E" = E t " intersecting at the same two-dimensional bifurcation sphere So- We call this sequence 
of slices E t (defined by the equation t = constant, with t' < t < t") a "tilted foliation" 0] . For a "tilted foliation" 
the spacetime domain M lying between E' and E" consists of two wedges A4+ and A4- located in the right (i?+) 
and left (R-) sectors of the Kruskal diagram. The region M. = A4 + U M.- is bounded by E' and E" and by a 
three-dimensional timelike boundary B that consists of two disconnected parts B + and For a general eternal 
black hole geometry the boundaries B + and B_ are located in M+ and respectively. We define S± =E± fl B±. 
The topology of the slices E is therefore Is x S 2 , where Iy, is a finite spacelike distance, while the topology of B± is 
It x i^., where It is a finite timelike distance, in Schwarzschild coordinates the line element of Eq.QJ can be written 



ds 2 = ~N 2 dt 2 



dy 2 



(y)dn 2 



(10) 



where y is the proper radial distance from the throat defined in Eq.Q. The coordinate y covers the entire range 
(— oo, +oo). N is the corresponding lapse function written in terms of the proper distance from the throat 



N 2 = 1 



2MG 
r(y) ' 



(11) 



We define the four-velocity vector with u M — —Nd^t, while the timelike unit vector normal to the hypersurface E t is 
u^. This is chosen to be future oriented in M. and normalized by the condition u ■ u = — 1. The lapse N is positive 
at E + , negative at E_ and vanishing at the bifurcation surface. The spacelike normal n M to the three-dimensional 
boundaries B±, is defined to be outward pointing at B±_- inward pointing at _B_, and normalized so that n-n = 1 with 
the further condition (u ■ n)\ B± — 0. Following Ref.[]J|, greek indices are used for tensors in M. while latin indices 
are used for tensors defined in either E or B±. The metric and extrinsic curvature of E as a surface embedded in M. 
are denoted by and i£y = —h^WkUj, respectively, while the metric and extrinsic curvature of the boundaries B± 
as surfaces embedded in M. are 7^ and 0y = — 7* : Vfcn J -. Covariant differentiation with respect to the metric and 
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hij is denoted by V and D, respectively. The induced metric and extrinsic curvature of the boundaries S± as surfaces 
embedded on £ are denoted by a a b and k a b = —a^Dkrib, respectively, (a, b = 2,3). Explicitly, the metric tensors for 
the different surfaces are 



The determinants of the metric tensors are related by 



(12) 



(13) 



The covariant form for the gravitational action generated by a tilted foliation with fixed three-dimensional boundary 
of M is 2 



S=7T- 



l - ( d i x v ^ g ~n-- [ d 3 xV^7(©- ©o)+ 

' K J M+ K JB+ 



1 

2k 



M- K 



d 3 Xy/^f(0- e )_ 



(15) 



where 5i denotes the four-dimensional scalar curvature, k = 8ttG. The integrations are taken over the coordinates x^ 
which have the same orientation as the canonical coordinates (t, y, 9, (f>) of the tilted foliation. The differing signs in 
M. + and M.- reflect the fact that the coordinates have different time orientations in M. + and M.-. The subtraction 
term 



S o(±) - ~ 



d 3 Xy/^f& 



(16) 



B± 



is the extrinsic curvature evaluated on the reference space, which in this case has been chosen to be flat. The effect 
of 5o(±) is to normalize the energy E to zero for the Schwarzschild wormhole with M = 0. In the next section we 
will carefully examine the contribution of the volume terms of the action. Here, we restrict the evaluation on the 
boundaries B± obtaining 



S = --( d 3 i v /z 7(e-6o) + - - / d 3 xy/^(e-@ ) 



The trace of the extrinsic curvature is 



6| ± =-V k n k = -- 



With the help of Eqs.© and JTT), we can write 

0|JV| 



-99 kj nj) = 



\N\r' 



■d y (\N\ 



dy 



± 



MG 

r 2 ^ 



in E H 



(17) 



(18) 



(19) 



where we have used 



dr 
'(hi 



= ±Wl- 



2MG 
r(y) 



in E±. 



(20) 
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The complete action should be considered with additional boundary terms of the type 



- / d?xVhK±. (14) 



K 

t' ,~E± 



Nevertheless, since we will look at the Euclideanized version of the action with periodic boundary conditions in the Euclidean time, the 
previous boundary term disappears. 
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Thus Eq.lJTSJ becomes 



Q \± = *\N\j2jfil 2r <3>)-3MG\, (21) 



while the subtraction term can be defined by 



o| ± = ^6| ± = T § (22) 



and 



(e - e„) ± = t 1 



|7V|r 2 (y) 



2r (y)-3MG-2 1 —^-\N\ 



(23) 



The proper "time" length in Ea. (|17|) is 



- f dt\N\ = -T\N(y±)\ in E±, (24) 

where y± is the location of the spatial boundary on B± and T = t" — t' . The substitutions i — > — zr and iV — * —iN E 
do not modify the behavior of the boundary action in the Euclidean context. Thus Ea. (|24ll becomes 

B ± =T ± 1 = N E {r ± )8* inS±, (25) 

where r± stands for r (y±) and where we have integrated over the period 8* = 8kMG obtained by imposing regularity 
condition on the horizon. Note that the inverse temperature at infinity is the same in both wedges. The difference 
between 8+ and /3_ essentially arises from the Tolman law involving the different location of the boundaries. Then 
the total action contribution is 

S B± = ± — T [2r - 3MG - 2rN E ±] | r=P± — ► I B± = T {8-kMt - 12irM 2 G - 0r/G) | r=r± . (26) 

Like in the case of the Hamiltonian discussed in Ref.^^j the tree level action for the Schwarzschild wormhole can be 
cast in the form 

- (1+ - /_) = (12itM 2 G - 8nMr + + P+r+/G) - (12ttM 2 G - 87rAfr_ + /3_r_/G) . (27) 

From Ea. (|27|l . we can compute the energy of the system at the fixed temperature f3± in each wedge. This is obtained 
by regarding Eq. l|27|l as a function of r and T rather than of r and M. We obtain 



01+ J_ 

00+ - G 

p _ ai- _ i 

Cj - — dp- — G 



r+ -r +v /l-2MG/r + 



-r-y/1- 2MG/r. 



(28) 



Note that each term tends asymptotically to the Arnowitt-Deser-Misner (ADM) mass Note also that if the 
boundary location is different in the respective wedge, implies that E + ^ E-. Nevertheless for the entropy we 
discover 



8+^-1+ = 0+E+ -I+ = 4ttM 2 G 
S- = f3- ^f- - I- = B-E- - J_ = 4ttM 2 G 



namely it is independent on the wedge. This is a consequence of the definition we have used for computing the 
entropy. The same result has been obtained by Martinez in Ref. [2(]] in the context of microcanonical approach to the 
entropy of an eternal black hole. Note that if one adopts the method used in Ref.[2(| for computing the entropy, it is 
not immediate to arrive at the result of Eq. 1)29(1 . By means of Eq. 1(29(1 , Eq. 1(27(1 can be explicitly rewritten as 

/3+E+ - (3-E- = 8ttMG {N e+ E + - N E -EJ) = 8* (M + - M_) . (30) 

Note that 8* is the periodically identified "asymptotic Euclidean time" which is independent on the considered wedge. 
We can recognize three different cases 
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1. M + > M_. The limit case is when r + — > oo (y + — > +00), while r_ = 2GAf (y_ = 0). This corresponds to 
Af + = A/^UjVj and M_ = 0. This situation is the standard one, where the black hole thermodynamics is 
investigated in the wedge M. + . Indeed 

2. M + < M_. The limit case is when r_ — ► 00 (y_ — » —00), while r + = 2GM (y + = 0). This corresponds 
to M_ = M^ujvi and M + — 0. This situation is "dual" to the previous in the sense that the black hole 
thermodynamics is investigated in the wedge M. _ . 

3. M_|_ = M_.The left boundary j/_ and the right boundary y + are symmetric with respect to the bifurcation 
surface Sq , which implies that r+ = r_ . 

Case 3 corresponds to a vanishing action. Nevertheless, it is well known that a vanishing action with a Euclidean 
metric describes flat spaceQ. Therefore, in this very particular case we have found an alternative way to flat space 
at the price of having a topology change. 



III. NONCONFORMAL NEGATIVE MODES 



The measure of stability in this approach is the eigenvalue spectrum of the nonconformal perturbative modes for the 
solution. Should there be a mode with a negative eigenvalue, then the action for this solution is a saddle-point in its 
phase space rather than a true minimum. Consequently, there ought to be a correspondence between the presence of 
such a negative mode and the local thermodynamic stability as governed by the heat capacity. Negative modes arise 
from the analysis of geometric fluctuations about classical Euclidean solutions of the Einstein field equations. However, 
the analysis to confirm their existence must be performed with care, since the gauge freedom of the Euclidean action 
will in general introduce a large number of non-physical negative modes associated with conformal deformations of 
the metric. For pure gravity, the contributions from the conformal and the nonconformal modes decouple if a suitable 
gauge is chosen. In the path integral approach [21|, the partition function Z is generally defined as a functional 
integral over all metrics with some fixed asymptotic behavior on some manifold A4, 

Z= f Z%]exp (-*%]). (31) 

JM 

This integral is formally defined by an analytic continuation to a Euclidean section of M. denoted by M. , where the 
right and left wedges of a Lorentzian eternal black hole (wormhole) are mapped into two complex sectors M.+ and 
M._ to become 

Z= f D[<7]exp (-/[<?]), (32) 

JM 

where the integral is performed over all positive definite metrics g. In the case of pure gravity, the Euclidean action 
comes from Eq. (|15|) whose expression is 

I = -^-l d 4 x^R+- [ d 3 x^(0-® Q ) 

ZK JM+ k JdM+ 



+ [ d 4 x^R+- [ d 3 x^(0- e ). (33) 

lK JM- K J'dM- 

This partition function may be approximated using saddle-point techniques, by Taylor expanding about the known 
stationary points of the Euclidean action - the solutions to the Einstein field equations 

= 0. (34) 

The expansions are performed by writing the metric g a b as 

9plv = 9^ + h^v (35) 

with h^y treated as a quantum field on the classical fixed background g M „ which vanishes on the boundary dM.±. 
Thus the Euclidean action becomes 



I[g] = I[g}+ + h[h]+ - I[g]- - h[h]- + 



(36) 
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where the linear term 1\ vanishes precisely because g^ is a classical solution, I2 is quadratic in the field hp U in each 
complex sector, and '• • • ' represents terms of higher than quadratic order. The Taylor expansion of the action in 
M. = A4 + U A4- leads to the following one loop approximation of the partition function 

logZ = -I [g] + log f D [h} + cxp (-I a [h} + ) + I [g}_ + log f D [h]_ exp (j a [h]_) . (37) 

JM+ J Ad- 

it is interesting to note that in this form the gravitational quantum fluctuations on A4 + LA/f_ J act as a sort of source 
field with respect to A4- (.M+). The quadratic contribution to the action is straightforward to evaluate, and may be 
written for arbitrary /i^„ in each wedge, in the form 

h[h]± = ^ j dtxy/gbgG^bg. (38) 

G is a physical gauge invariant spin-2 operator which acts on the transverse and trace-free part of h^ v and takes the 
simple form 

G^ypa — 9 tip9vo^j 'qV ^R^pua • (39) 

In general, G will have some finite number - generally zero or one - of negative eigenvalues, which correspond to the 
nonconformal negative modes of the solution. The eigenvalues of G arc determined by all solutions to the elliptic 
equation 

\n)h% (40) 

where the eigenfunctions hp V are real, regular, symmetric, transverse, trace-free, and normalizable tensors. Clearly, 
should one of the eigenvalues of G be negative, then the product of all of the eigenvalues would also be negative. 
The contribution to logZ from fluctuations about the classical solution would then contain an imaginary component, 
leading to an instability in the ensemble similar to the type proposed by Gross, Perry and Yaffe Naively, then, the 
one-loop term may be written as \ log Det(/i~ 2 G) where \x is a regularization mass, and the determinant is formally 
defined as the product of the eigenvalues of G. However, due to the diffcomorphism gauge freedom of the action, A 
will in general have a large number of zero eigenvalues, and so this procedure as stated is ill-defined. The remedy 
is to add a gauge fixing term B - such that the operator A + B has no zero eigenvalues - and an associated ghost 
contribution C, to obtain 

logZ = -I[g] - i log Det(^ 2 {A + B}) + log Det{^ 2 C). (41) 

Such terms may be dealt with by means of generalized zeta functions, as considered by Gibbons, Hawking, and 
Perry , and extended to include a A term by Hawking In order that this be possible, the terms must be 

expressed as sums of operators, each with only a finite number of negative eigenvalues. This may be achieved by 
writing A + B as —F + G, where F is a scalar operator acting on the trace of h^. The ghost term C is a spin-1 
operator acting on divergence-free vectors, and An observation of primary significance is that a gauge may be chosen 
in which the operators F and C have no negative eigenvalues. If the background metric g a b is flat then, in addition, 
G will be positive-definite, but for a non-flat background this is not the case. 



IV. VARIATIONAL APPROACH TO THE NEGATIVE MODE 



Following the treatment given in 4], and subsequently in |6( for the Schwarzschild case, it is clear that only 
spherically symmetric and r-independent solutions of Eq. I|40l) need be considered as candidate nonconformal negative 
modes. In static and spherically symmetric backgrounds, modes of higher multipole moment will necessarily have 
greater eigenvalues. With this assumption, it is then straightforward to write down a construction for such solutions 
to G valid in a four-dimensional Euclidean wormhole background of the form (|10|l . Since G acts only on symmetric 
transverse and trace-free tensors, then clearly the constructed solutions must exhibit all of these properties. If the 
mode h£ v is written in the manifestly trace-free and symmetric form 

(Kf = diag (H±{y),Ht{y)-\ (H±(y) + ff±(j/)) (iJ±(y) + H±{y))\ , (42) 
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then the final property, V a /i*„ = 0, is guaranteed if it is further assumed that Hq(u) and H^(y) are related through 
the first order equation 



3 , N v r , N v + 3 , 



(43) 



In Eas. i|42[l and l|43() . we have explicitly written the location of the perturbation on each patch covering one universe. 
The two patches join at the throat of the wormhole and the transversality condition can be cast into the form 



dHf{r) 3 



1 dN 



dr 



-Hf{r) + -—H?(r) 
r N dr 



dr 

dy 



Ndr H ° [r) r olJ 



dr 



(44) 



where it is manifest the independence on the wedge location. This is a consequence of the discrete symmetry (JJJ 
between the asymptotically flat wedges. When we apply this symmetry to the the eigenvalue equation (|45|l with the 
help of the following substitution 



N 



-N 

~T .1/ . 



we find that Ea. l|45[) is invariant in form only on the coordinate r. With the ansatz (|42|) and 144(1. the eigenvalue 
equation (14011 reduces to a linear second order ordinary differential equation for the component H f' (y) and eigenvalue 
fi 2 . To further proceed, we insert Eqs.10 and into Eq. l|4U|l . Since we have chosen y and the Euclidean lapse Ne 
to be positive in 7W + and negative in A4-, then Ea. (|40|l is separated in two pieces 3 



d 2 Hi (y) , ,± , dHi (y) , ± f , tt / s 

+■ b (y) — i — + a (y) Hi (y) 



dy 2 



dy 



where 



b + (y) = Wi 



in M+ and 



b~ (y) 



2 MG (4 r 2 (y) - 23 MGr (y) + 27 M 2 G 2 ) 
' r (y) r 2 (y) (r (y) - 2 MG) (r (y) - 3 MG) 

(4 r 2 (y) - 23 MGr (y) + 27 M 2 G 2 ) 
' r 2 (y) (r (y) - 3 MG) y/T^2MG/r~(y) 



-^Hx (y) , 



a + (y) = 



8MG 



(r (y)-3MG)r 2 (y) 



a (y) 



8MG 



(r (y)-3MG)r 2 (y) 



(45) 



(46) 



(47) 



in M.-. After having cast Ea. (|45|l in the standard Sturm-Liouville form, we apply a variational procedure to establish 
the existence of an instability by means of trial functions. The standard Sturm-Liouville form of Ea. H45|) is 



dy V d y 

where the integrating factor (y) is defined by 
P ± (y) = exp / b ± (y) dy = exp ± 



P ± (V) -rH, (y) + q ± (y) H x (y) + fi 2 p ± (y) H x (y) = 0, 



b ± (r(y)) 



y/l-2MG/r (y) 
To compute the integrating factor we used Eq.©. The coefficient 

q ± (y)=P ± (y)a ± (y) 

and the boundary conditions 



zdr = 



(r(y)f 2 (r(y)-2MGf 2 
(r(y)-3MG) 2 



V (y) — j Hi (y) 



dy 



and 



v=o 



p yy) — — h 1 ( y ) 

dy „=_, 



(48) 



(49) 



(50) 



(51) 



3 The analogue of Eg. 1451 has been studie d i n the context of the s tabi lity of th e S chwarzschild solutio n fo r the first time by Regge 
and Wheeler^, VishveshwaraQ , Zerilli^, Press and Teukolsky, ^j^, Stewart |27^ and ChandrasekharQ. They concluded that the 
Schwarzschild solution is classically stable. 
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lead us to write the Sturm-Liouville problem in the following functional form 

2 



[Hl {V) ; A] = D[ Hl (vhX] = M± (A) = J I±d y P Hy)H!(y) ' (52) 

where I + = [0, +00) and i_ = (— 00, 0]. The choice of a trial function is suggested by the asymptotic behavior of in 
Ea.(|45p. Indeed, when y — ► ±00, (y) ~ exp (=fAj/) ~ exp (— Ar). The denominator in Eq.JS^J) compensates the 
lack of normalization in the trial function H± (y). If A is defined as a variational parameter, then the Rayleigh-Ritz 
method leads to the following result 4 

M± (A min ) = -.1374637652 for A min = .4471243737. (53) 

However the denominator of Ea. (|52f) vanishes for A = 1.360771766 and for A > A, the 11 normalization" changes sign. 
The lack of positivity for every A is related to the appearance of a spurious singularity in r = 3MG 5 . Therefore, it is 
necessary to choose a better form for the trial function (y) which partially removes the singularity in r = 3MG. 
Such a form is given by H l (y) ~ Hi (r) = exp (— Ar) (r — 3MG). With this choice, the "normalization" is strictly 
positive for every A. The minimum of fi± (A) is reached for A min = .9223711858 and n\ (A min ) = —.3827139389. It is 
interesting to note that 

f 2 (A min ) - /x+ (Anxin) + fil (A min ) = -0.765 44 - M 2 (A min ) / (2MG) 2 = -0.191 36/ (MG) 2 . (54) 
This is the same value obtained by GPY for a single Euclidean black hole. 



V. SUMMARY AND CONCLUSIONS 



In this paper we have analyzed, the one loop contribution coming from quantum fluctuations around a background 
metric describing a Schwarzschild wormhole which has boundary terms in both wedges of the Penrose diagram. 
Usually, this kind of analysis is performed only in the R+ wedge of diagram l|T|l. However, due to the symmetry 
property of transformation (0, a contribution from region i?_ is not unexpected. Moreover, the possibility of a fine 
tuning between boundaries in the opposite regions gives the opportunity to vanish the classical gravitational action. 
The vanishing of the classical action leads to the interpretation that we are dealing with flat space. On the other 
hand the presence in each wedge of a term proportional to the ADM mass is a signal that a non flat configuration 
has been taken under examination. An interpretation of this puzzling situation can be given in terms of black hole 
pair creation mediated by a wormhole, with the pair elements residing in the different universes. This proposal is 
not completely new and it has been investigated in Refs.fl6j| in a Hamiltonian formulation, where a foliation of the 
manifold M. is crucial. On the other hand, the action formalism has the advantage of being fully covariant. However, 
the vanishing of the classical action is not sufficient to establish if spacetime will produce a pair or will persist in its 
flat configuration. It is, therefore necessary to discover if the one loop approximation will produce an instability. The 
variational approach we have used for the Sturm-Liouville problem in the TT sector has shown an unstable mode in 
the s-wave approximation in both wedges. This is the signal that flat space is not only unstable when the space is 
furnished with a temperature T = 1/ (8irMG), but is even unstable with respect to Schwinger pair creation. One 
possible conclusion is that flat space can be no more considered as the general accepted vacuum of quantum gravity 
and the instability with respect to Schwinger pair suggests to search in this direction instead of rejecting the result 
with the further motivation that this instability can disappear if the bounding box is restricted enough. A partial 
success in this direction has been obtained by Allen 6] in terms of hot quantum gravity from one side. On the other 
side, the choice of a non-trivial vacuum formed by a large number of black hole pair generated in the same way we 
have investigated in this paper has led to a possible stable picture of space time foam, which could replace flat space 
as a candidate for the quantum gravity vacuum p9j. 
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4 Details of computation can be found in the Appendix. 

5 See details in Appendix lAl 
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APPENDIX A: SEARCHING FOR EIGENVALUES WITH THE RAYLEIGH-RITZ METHOD 



In this Appendix, we will explicitly solve the eigenvalue problem of Eq.{a2} with the Raylcigh-Ritz method. We 
here report two principal choices of trial wave functions: 



1. iff (y) = cxp (^fXy) — ► exp (— Xr). For this choice Ea. H52(l becomes, 



where 



p 2 (A) = A — 8MG— 



r + °° dre~ 2Xr r ( r - 2MG ) AIV , 

hue'-" (i — 3MG) _ ^2 A W 



I2MG ^ re 



-2\ r r 5 (r-2MG) 
(r-3MG) 2 



B(X) 



A (A) = 8MG 



2MG 



dre _ 2Ar r3(r-2MG) = ^ ^3 ^ ^ pHp -1) 



(r - 3MGY 



(2p-3) a 



(Al) 



and 



256 (MGY 



It ( A ) - /! ( A 



S(A) 



2MG 



dre ~ 2Xr V^GY = 128 (MG)5 f °° dpcxp ^ 



p 5 [p - 1) 

(2p-3) 2 



(A2) 



= 128 (MG) 5 if (a) - if (A 
We have defined p — rj (2MG), A = X2MG and introduced the following function 



/" A 



d n Il X 



dpexv[-2Xp) =(-!)". 



dA r ' 



4n f^j can be easily integrated with the formula |30l| 

x) / m — 1 l\\ ( \in — fe— 1 / \ vn 1 

^ = ^P (-H V } V, jr- - V TTi exp [fiv) Ei [- ( u + 0) u] 

(x + P) fe( m_1 ) ! + (m-1)! MJ 

with to > 2, |arg (it + /?) | < 7r and Re v > 0. Ei (x) is the exponential function. We get A (A) = 

(-18 A e A - 108 A 3 Ei (A) + 72 A 2 Ei (A) + 81 A 2 e A - 27 A 3 e A + 27 A 4 Ei (A) - 2 e A ) e~ 3 A 

8A 2 

and B (A) = 

(24 e A - 243 A 5 e A + 243 A 6 Ei ( A) - 648 A 5 Ei (A) + 194 A 3 e A + 72 A e A + 126 A 2 e A + 405 A 4 e A ) e~ 3 A 

128A 5 ' 

B (A) vanishes for A = 1.360771766 and p 2 (A) gets its minimum for A m j 
-.1374637652 



(A3a) 



(A4) 



(A5) 



(A6) 



(A7) 



.4471243737 leading to p 2 (A min ) 



2. With the help of Ea. (|X5l) and choosing Hf (y) = exp (TXy) (r (y) - 3MG) -> cxp (-Xr) (r - 3MG) we get a 
more complicated form in N [Hi (y) ; A] of Eq. l|52|l . while D [Hi (y) ; A] is simpler. Indeed 



r+oc 

N [Hi (y) ; A] = - / dr 

J2MG 



r 5 (r-2MG) ( d r _ 2Xr 



2MG (r — 3 MG) \dr 



)2 r+oo 
+ 8MG / dre 
J 2MG 



_ 2Xr r*(r-2MG) 
(r - 3MG) 



(2MG) 7 f + ™ d pHp-1) ( d 
(2MG) 2 J i P {2p-3) 2 \dp 



e- 2 ~ xp (2p-3) 



8MG 



(2MG) 5 
MG 



dpe 



(2p - 3) 
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(2MG) 



(243 X 6 Ei (A) - 239 A 5 e A + 42 e A + 225 A 4 e A - 432 X 5 Ei (A) + 86 A 2 e A + 102 A 3 e A + 78 A 



e A )e- 3A 



16A 5 

(A8) 

and 



D (y) ; A] = / dr [e" 2A ' V 5 (r - 2 MG)] = (2MG) 7 / dp exp ( -2Xp ) p 5 (p - 1) 

J2MG Jl ^ ' 

e" 2A (45 + 2 A 5 + 10 A 4 + 30 A 3 + 60 A 2 + 75 A) , 
= (2A/G) 7 ^ ^ i . (A9) 

Note the positiveness of D [Hi (y) ; A] for every A, ensuring a correct normalization. The minimum of /i 2 (A) is 
reached in this case for A min = .9223711858 and p 2 (A min ) = -.3827139389. 
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